Transmutation from Skyrmions to Half-Solitons 
Driven by the Nonlinear Optical Spin-Hall Effect 
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We show that the polarization domains generated in the linear optical spin-Hall effect by the ana- 
log of spin-orbit interaction for exciton-polaritons are associated with the formation of a Skyrmion 
lattice. In the non-linear regime, the spin anisotropy of the polariton-polariton interactions results 
in a spatial compression of the spin domains and in an abrupt transformation of the Skyrmions 
into half-solitons, associated with both the focalization of the spin currents and the emergence of a 
strongly anisotropic emission pattern. 
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Introduction. The ultimate goal of spintronics is to 
replace charge currents by pure spin currents. The spin- 
Hall effect is one of the key concepts in this research field, 
since it allows the formation of a spin current perpen- 
dicular to the electric current without a magnetic field. 
Originally proposed by Dyakonov and Perel' in 1971 [l| 
and rediscovered by Hirsch in 1999 [2], it attracts a lot 
of attention now 0, 0] • Basically, this effect originates 
from the electron's spin-orbit interaction that makes its 
scattering on impurities spin-dependent. One of the most 
serious obstacles to the realization of spintronic devices 
is the dramatic role played by the processes of spin re- 
laxation, inducing information losses. The optical coun- 
terpart of spintronics, namely spin-optronics, has been 
proposed as a valuable alternative, since the characteris- 
tic decoherence times of optical excitations are orders of 
magnitude longer than those of electrons and holes [5[. 

Exciton-polaritons (polaritons), being a mixture of 
quantum- well excitons and cavity photons, appear as the 
optimal candidates for the realization of optoelectronic 
devices showing a polarization-selective behavior. The 
large photonic fraction of polaritons allows an efficient 
resonant optical excitation with a laser pump at any en- 
ergy, momentum and density, with a well defined spin to 
form macroscopically occupied and coherent states hav- 
ing a controllable spatial extension. Moreover, the po- 
lariton mass is five orders of magnitude smaller than that 
of free electron, allowing their propagation at large ve- 
locities over hundreds of microns |6[ preserving the spin 
coherence 0- On the other side, their excitonic part al- 
lows them to interact efficiently with each others. These 
interacting photonic particles are therefore demonstrat- 
ing an outstanding non-linear optical response [!, Q . 

Polaritons have two spin projections on the structure 
growth axis of the sample [lflj. An analog of the spin- 
orbit interaction is present in microcavities, provided by 
the so-called TE-TM splitting [ll|. This feature gives 



birth to the optical spin-Hall effect (OSHE) [12| which 
became a popular research topic [13l-ll6j. The OSHE 



however remains a linear effect, which does not capital- 
ize on the exceptional non-linear optical response of the 
polar iton system [T^j]. Furthermore, the created spin cur- 
rents are by nature relatively weak, and demonstrate an 
enormous angular aperture of almost 45 degrees. 

Polariton-polariton interactions are strongly spin 
anisotropic [17|. which imposes a spin-dependent nonlin- 
ear response |l8l~[2o|. resulting in the formation of spin 
polarized topological excitations namely half-vortices 
[2l|,l22[ and half-solitons [23|,l24|. These objects have cir- 
cularly polarized dips on a linearly polarized background 
[25I ]. the density in the center being one half of that at 
infinity. They differ from the topological excitations in 
spin-isotropic condensates, which are Skyrmions or core- 
less vortices [26], appearing as domains of circular po- 
larization inversion, while the total density remains con- 
stant. It was recently understood [27] that half-integer 
topological defects in polariton condensates behave as 
stable magnetic charges accelerated by effective magnetic 
field. They can therefore be foreseen as ultrafast spinop- 
tronics information mediators [28| . 

In this letter, we describe the spectacular transition 
arising between the linear and nonlinear regimes in the 
OSHE. In the linear regime, we demonstrate that the 
spin domains generated by the interplay between the 
particle flux and the spin-orbit like interaction are car- 
rying topological charges, similar to Skyrmions forming 
a lattice. In the nonlinear regime one observes a dra- 
matic change of the polarization pattern, emission profile 
and nature of the topological excitations. The spin do- 
mains are found to be spatially compressed by the spin- 
anisotropic polariton-polariton interaction. The phase 
transition occurs when the healing length of the nonlin- 
ear fluid becomes smaller than the size of the polarization 
domains. The Skyrmions created by the spin-orbit inter- 
action abruptly transform into half-solitons, stabilized in 
turn by the interactions. Finally, the total emission pro- 
file switches from isotropic to strongly anisotropic, ac- 
companied by a very strong focalization of the spin cur- 
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FIG. 1: (a) Pumping scheme showing the lower polariton 
branch (LPB) and its excitation on a ring by the pump spot. 
Inset: cylindrically symmetric distribution of the total inten- 
sity under the analytical description, (b) Degree of circular 
polarization p c showing the antisymmetric circular polariza- 
tion domains, (c) Zoom on the phase of the cr+ component 
characterized by a fork-like dislocation at the position of the 
density minimum which is the typical signature of a topolog- 
ical defect in this component. No dislocation occurs at the 
<j- minima (red regions of panel in (b)). (d) Pseudospin tex- 
ture zooming on a Skyrmion and showing the winding of the 
in-plane pseudospin around its core (black arrows). 



rents. 

Linear OSHE. The original proposal for the OSHE [12 ] 
involved the resonant injection of exciton polaritons with 
a well-defined momentum and linear polarization within 
the disorder landscape of a microcavity. The resulting 
coherent Rayleigh scattering allowed a homogeneous re- 
distribution of polaritons on an elastic ring in /c-space. 
However, in modern microcavities the disorder tends to 
be strongly reduced and would anyway be screened by 
the particles interactions in the nonlinear regime. An al- 
ternative configuration for the observation of OSHE is to 
use a blue detuned laser with a small spot [l4| giving a 
sufficient extension in momentum space to excite directly 
the polariton dispersion on a ring [see FigQJa)]. The 
k-dependent precession of the polariton pseudospin S 
around the effective magnetic field Hlt associated with 
the TE-TM splitting results in the appearance of alter- 
nating circularly polarized domains in the four quarters 
of the plane, as shown in FigJTJb). 

For simplicity let us consider the ideal situation of 
particles continuously injected at r = 0, linearly polar- 
ized along the x direction and well defined energy. This 
boundary condition excites a specific linear combination 
of TE and TM eigenmodes for each polar angle. The 
resulting polariton wave function, stationary solution of 
the spin-dependent Schroedinger equation with external 



pumping term written here on the (x, y) polarization ba- 
sis reads 



e iktT cos 2 p + e iklT sm 2 (p) 



sin ip cos p + e 



(i) 

' sin cp cos ip) 
(2) 



Here ro = hkr/m* with r and m* are the polariton life- 
time and effective mass respectively and k = (ki + kt)/2. 
kij = ^2m*Ei )t /h where E^ t is the energy of the TM 
and TE mode respectively. The distribution of the total 
emission shown on the inset of the FigOJa) is assumed 
here to be fully isotropic. FigQJb) shows the degree of 
circular polarization p c = ^ (ip x ip y ) / (n x + n y) of the po- 
lariton emission. It demonstrates the typical periodic po- 
larization oscillations with maxima lying along diagonal 
directions, perfectly reproducing the polarization pattern 
of the original OSHE. The positions of the maxima are 
defined by 
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where p G Z and the indices ± refer to the <j± com- 
ponent respectively given that i/;± = {$ x ± iip y )/y/2. A 
fundamental feature is that the cores of the spin domains 
are characterized by phase singularities in the circular 
polarized components shown in FigOJc). They appear 
because the effective field projections on the x axis are 
exactly opposite on the two sides of the diagonal direc- 
tion ip = £>7r/4, inducing a tt phase difference after a 
half-period rotation of the pseudospin. This is the time- 
reversal of the half- vortex unwinding scheme [29| , which 
in our case produces not a half-vortex, but a Skyrmion 
excitation, characterized by a homogeneous total density 
distribution [see inset of FigJTJa)]. The latter appears 
because the spin- anisotropic interactions, known to sta- 
bilize the half- vortices, are absent yet, and the system is 
governed instead by the spin-is otropic effective magnetic 
energy. The pseudospin texture zoomed on a Skyrmion 
is shown in FigJTJd). 

Nonlinear OSHE. When the density increases, the 
quantum fluid becomes non-linear, and can support topo- 
logical excitations driven by the interparticle interac- 
tions, their size being defined by the healing length £. 
Due to the spin anisotropy of the polariton interactions, 
Skyrmions cannot appear any more and are replaced by 
half-solitons. Qualitatively, one can expect the transition 
to occur when the healing length becomes smaller than 
the Skyrmion size (of the order of r^). Once transformed, 
the resulting defect starts to feel the effective magnetic 
force provided by the TE-TM splitting, which now acts 
as a weak perturbation. 

The complete simulation of the problem requires nu- 
merical solutions of non-linear wave equations. We model 
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the system of resonantly injected cavity polar itons with a 
system of equations for the coupled photonic <fr(r, t) and 
excitonic fields x( Y ->^)-> each of which have two compo- 
nents written here on the circular polarization basis: 
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In the above expressions Vr = 5 meV denotes half of 
the Rabi splitting coupling excitonic and photonic fields, 
r x = 400 ps and = 10 ps are the lifetimes of excitons 
and cavity photons respectively, P(r,t) is a 2 fim large 
Gaussian, quasi-resonant and linearly polarized (along 
the x-axis) optical pump having a frequency uop blue de- 
tuned by +1 meV from the lower polariton branch. m x = 
0.4rao, rriff) = 5 x 10 -5 rao are the effective masses of the 
excitons and cavity photons respectively (mo is the elec- 
tron mass). The constant j3 = H 2 /4(l/m% M - l/m% E ), 
where rn^ E = and m^ M = O.95m0 are the masses 
of the photonic TE and TM modes, is the strength of 
the photonic TE-TM splitting. The corresponding terms 
give rise to the in-plane effective magnetic field H^T(k). 
The constants ql\ and ot2 = — 0.2ai describe the spin- 
anisotropic exciton- exciton interactions. 

The results are presented in Figj2j Panels (a) and (e) 
show the circular polarization degree in the linear and 
nonlinear regimes, correspondingly. In the linear regime 
one can notice that the spin pattern is not symmetric, 
with the spin domain being slightly compressed along 
the y-axis. Indeed, in a realistic simulation (and experi- 
ment) the use of a small size, monochromatic laser tends 
to excite TE and TM modes with two different wave 
vectors instead of two different energies as assumed in 
the original OSHE and to derive Eqs.(l,2). This induces 
two different velocities for the TE and TM propagating 
waves which would result in two different effective decay 
lengths. This is slightly shifting in (c) the directions of 
maximum circular polarization away from the diagonal 
direction and slightly breaks the homogeneous distribu- 
tion of the total intensity. As soon as the polariton den- 
sity increases so that the healing length £ ~ h/ ' ^Ja\n x m§ 
becomes smaller than the size of the Skyrmions deter- 
mined by , the nonlinear effects start to dominate over 
the TE-TM splitting, and the interplay between the two 
leads to the transmutation (conversion) of Skyrmions into 
half-solitons. The latter are characterized by phase dis- 
locations clearly visible close to the x-axis in (h). They 
also appear in (f) as a domain walls between x and y 
polarizations [23| and in (g) as sharp minima in the to- 
tal density distribution. The angular dependence of the 



emission is therefore switching from smooth to strongly 
anisotropic. Interestingly one can see in (e) an oscillation 
of the circular polarization along the half soliton trajec- 
tory. Indeed, the decay of the total density far from the 
pump spot makes the half-solitons unstable against the 
effective field Hlt as soon as it exceeds the critical field 
H c = (oi\ —Oi2)nj^ [I?}. The two half-solitons of opposite 
circularity therefore get coupled, leading to spatial oscil- 
lations of the circular polarization degree. However, due 
to the finite life time of polaritons the integrated emission 
the solitons axes remains strongly circularly polarized, as 
visible in FigEJb). 

In addition, we have extracted in Fig|3]^a) the angular 
extension Aa(cp) of the topological defects (solid/purple 
line) and their local contribution to the total density 
n re i = 1 — nj/(nj + n^) (dashed/blue line) as a function 
of the pump spot density no- The abrupt evolution at a 
critical density of both quantities demonstrate that we 
assist to a phase transition. Indeed the angular aperture 
of the topological defects switches from 40 degrees to less 
than 5 degrees and the relative density from to almost 
0.5 where it remains locked. 

The compression of the spin domains can be described 
in the framework of the pseudospin dynamics solely [To| , 
which cannot address however the formation and dy- 
namics of topological defects. The pseudospin preces- 
sion equation reads d t S = H x S/H. The total effective 
magnetic field H represents the sum of the field provided 
by TE-TM splitting Hlt lying in plane of the cavity 
and the effective field provided by the polariton-polariton 
interactions, responsible e.g. for the self-induced Lar- 
mor precession [30[ (intrinsic Zeeman splitting). The 
latter is directed perpendicular to the plane of the mi- 
crocavity (z direction) and describes the splitting be- 
tween the circular polarized polariton states appearing 
due to the spin anisotropy of the interactions: = 
— (ql\ — 012) (n+ — ri-) /2u z . The resulting system of non- 
linear differential equations cannot be solved analytically 
in the general case. 

If the pseudospin is initially aligned along the x-axis, 
important conclusions can be made studying the dynam- 
ics at short times. Indeed, if one considers the dynamics 
of the S y pseudospin projection (diagonal polarizations), 
it turns out that the contributions of the TE-TM and 
nonlinear fields to this projection are able to cancel each 
other if the propagation direction angle cp verifies the con- 
dition Ho cos(20) = —an. This compensation determines 
the direction of the maximal oscillations of the polariza- 
tion degree: instead of rotating around the in-plane field, 
the pseudospin turns further and further towards the z 
axis, because its S y projection remains zero. In the linear 
case, this direction corresponds to the TE-TM field com- 
pletely orthogonal to the initial pseudospin <p = pir/A, 
while in the non-linear case the maximum is shifted to- 
wards the y axis, along which the field is antiparallel to 
the injected pseudospin, until the size of the defects be- 




FIG. 2: (a-d) Linear and (e-h) nonlinear regimes. (a,e) Degree of circular polarization showing the compression of the domains 
and the half-soliton formation (e) while the density is increased. (b,f) Degree of linear polarization pi = (n x — n y )/(n x + n y ). 
(b) the linear polarization oscillates in the oblique direction, (f) The half-solitons are clearly evidenced as domain walls between 
x (red regions) and y (blue regions) polarizations. (c,g) Total emission, normalized to compensate for the global decay, goes 
from almost homogeneous (c) towards the appearance of visible minima at the position of the half-solitons. (d,h) Interference 
patterns: In the linear regime (d) only the cr+ component is plotted: typical forklike dislocations at the position of the density 
minima evidencing topological defects appearing in only one component, (c) The homogeneity of the total emission means that 
these defects are Skyrmions. (h) Phase of the total emission revealing an elongated phase slip along the total density minima 
evidencing further the half-solitons. 



comes smaller than r^, when the phase transition occurs. 

The focalization of the spin currents is illustrated in 
Figj3]^b) showing the analytical (solid lines) and numeri- 
cal (dashed lines) space-integrated degree of circular po- 
larization p c ((f) = J dr(n+ — n_)/(n + + n_). A sub- 
stantial compression of the spin domains is obtained. In 
the linear regime (black lines) we see that the spin cur- 
rents are almost unusable while in the nonlinear regime 
(red lines) the focusing is promising for a real spin cur- 
rent generator. The focalization however is much smaller 
in the analytical case than in the full Gross-Pitaevskii 
model [Eqs.flnEJ)], shown in Fig|3]^b), and which takes 
into account the formation and dynamics of topological 
defects. 

Conclusions. In summary, we have found that the lin- 
ear optical spin-Hall effect is associated with the forma- 
tion of a Skyrmion lattice. In the nonlinear regime a 
transmutation towards pairs of half-solitons occurs. The 
latter are associated with the breaking of the cylindrical 
symmetry of the total density and the strong focusing 
the spin currents in the real space. 
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FIG. 3: (a) Phase transition: angular extension of the topo- 
logical defects Aa (solid/purple line and left scale) and their 
relative contribution to the total density n re i (dashed/blue 
line and right scale), (b) Focalization of the spin currents: 
Analytical (solid lines) vs numerical results (dashed lines). 
Black (red) curves: spatially integrated degree of circular po- 
larization ~p c in the linear (nonlinear) regime. 
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